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I. INTRODUCTION 

PLANETRACK is a program in the Haystack Univac 490 pointing systeni whose 
output is the celestial coordinates of a selected planet at a given time. The progran; 
computes the coordinates by 3rd or 4th difference interpolation in the tables of the 
apparent right ascension, declination, and radius vector of the planets published 
in The American Ephemeris. The rates of change of the three quantities are computed 
by numerical differentiation. 
Il. INPUTS TO PROGRAM 

A. Inputs Furnished by Core Memory 

The program uses the registers listed in Table I for input information. The year 
and day are used to select the appropriate entries from the ephemeris. The coordin- 
ates of the planet are interpolated for the time which is the sum of the times in CEL- 
TIME and DELTATEE. 

B. Inputs Furnished by Magnetic Tape 

The tape format of the Ephemerides for the planets, which has been edited and 
recorded on tape by a 7094 program described in Reference (1), is shown in Fig. 1. 
Each block of 320 words covers a period of 32 days. The ephemeris for each day 
requires 10 words. The program does not make use of the semi-diameter, horizontal 
parallax, or the year-month-day words in the block. The tape must be on Unit 1 
(normally Servo B), of the Univac 490. 

C. Inputs Furnished by Typewriter 

During the operation of the initialization section of PLANETRACK an identification 
number representing the selected planet must be entered via the typewriter. 

II. PROGRAM OUTPUTS 

A. Outputs Left in Core Memory 

The apparent celestial coordinates of the planet, corresponding to the time in the 
register CELTIME, together with their numerical derivatives, are stored in the 


appropriate registers, as illustrated in Table II. 


Register Contents and Scaling 


W(CELTIME) Days B28 
W(DELTATEE) Days B28 
U(YEARMONTH) Year B15 (4 decimal digits) 
L(YEARMONTH) Month BO 
L(DAY) Day of Year BO 
TABLE | 


Core Memory Inputs to PLANETRACK Program 


Note: Notation "B28 means that the 


binary point is to the right of bit 28. 


(Planet number) 


4 8 4 Q) 
Right Ascension B26 Radians 


First Difference Right Ascension B26 Radians 


Declination B26 Radians 
First Difference Declination B26 Radians 
Radius Vector B24 Astronomical Units 


First Difference Radius Vector Astronomical Units 


Semi-Diameter B26 Radians 


Horizontal Parallax B26 Radians 


Fig. 1 Magnetic Tape Format for Planet Ephemeris 


Register Contents and Scaling 


W(RA) Apparent Right Ascension B27 in revolutions 
W(DEC) Apparent Declination B27 in revolutions 
W(RADIUS) One's Complement Distance from center of 


Earth to Surface of Planet B24 in Astro- 
nomical units 


W(RADOT) Numerical Derivative of Right Ascension 
B37 radians/sec. 

W(DEDOT) Numerical Derivative of Declination B37 
radians/sec. 

W(RADIUSDOT) Numerical Derivative of Radius Vector B24 


nautical miles/sec. 


TABLE II 


Core Memory Outputs of PLANETRACK Program 


B. Typewriter Outputs 
When the PLANETRACK program is initially called, it will use the INTERCOM 
subroutine to request from the operator the name of the planet whose coordinates are 
required, The typewriter will print the names of the planets, (except Earth and Pluto) 
each followed by the identification number in parenthesis. 
During initialization the PLANETRACK program will type, using the INTERCOM 
Subroutine, certain information concerning the selected planet. ltems printed are: 
lL. Julian Day corresponding to values in YEARMONTH, DAY, and CELTIME 
(7 digits). 
2. Apparent Right Ascension of planet (Hours, Minutes, Seconds to hund- 
redths). 
Apparent Declination of planet (Degrees, Minutes, Seconds to hundredths). 
Day of Year corresponding to values in YEARMONTH and DAY (up to 
3 digits). 
5. Universal Time for which the coordinates are interpolated. It is the time 
in CELTIME at initialization (Hours, Minutes, Seconds to hundredths). 
6. The distance between the centers of the Earth and planet (Astronomical 
Units to hundred millionths). 
7. The name of the planet. 

lf, after a search through 8 files on Unit 1 (or the finding of an end of tape mark), 
the planet ephemerides are not found, the program will type "(W) EPHEMERIS FOR 
X/Y NOT AMONG FIRST 8 FILES" and will exit to the master control program error 
return. ‘W"' is the number specifying the selected planet, "X" and "Y" indicate the 
current month and year, respectively. 

If there is a tape servo malfunction during search, the program will type "IIIC 
STATUS Sl ZZ" and will exit to the master control program error return. The two 
digit octal number "ZZ" comes from the tape status word and indicates the error con- 
dition. 

IV. ASTRONOMICAL SIGNIFICANCE 

The apparent right ascension and declination of the planet in the emphemerides 
are referred to the true equinox and equator of date and are corrected for planetary 
aberration. They are geocentric apparent quantities; the parallax correction is made 


by the coordinate conversion program (COCON) in the Haystack system. The values 


in the tables are computed for Ephemeris Time as argument. The equation 
ET = UT+AT 


is used to convert from Universal Time to Ephemeris Time. The constant AT, 
represented by the number in the register DELTATEE is approximately 35 seconds. 
The radius vector is the actual geometric distance in astronomical units between 
the centers of the Earth and planet at the stated time. 
The number 934.91 is used to convert astronomical units per day to nautical miles 


per second in the computation of the number in RADIUSDOT. 


499. 005 x 2.997925 x 10° 


934.91 = 1. 852 x 86400 


499.005 = number of seconds per astronomical unit 
2.997929 & 10° = velocity of light in kilometers per second 
1.852 = kilometers per nautical mile 


86400 = seconds per day 


A comprehensive explanation of the ephemerides is given in Reference (2). 
V. PROGRAM DETAILS 

The PLANETRACK program is a subroutine of the Haystack Univac 490 pointing 
program. ‘The initialization section begins at PLANETINIT, the working section at 
PLANETCONT. ‘There are several closed subroutines within PLANETRACK. ‘These 
are: DAYFIND, STATUSCK, INTERPOL, LEFRNDOFF, and ROUNDOFF. A listing 
of the program is given in Appendix I. 

A. Initialization 

The program, upon initialization, stores an RIL instruction in the tape channel 
internal interrupt register and an RJP STATUSCK instruction in the external interrupt 
register. ‘The area in core into which the tape data will be read is cleared, together 
with some additional registers. ‘This is done to make diagnosis easier in case of 
malfunction. 

The planet ephemeris entries are serially numbered by the program which gene- 


rates the magnetic tape (Reference 1). These serial numbers are called ‘Conventional 


Day Numbers" (CDN). An arbitrary decision was made to produce and use tapes such 
that the CDN for 25 April 1963 is zero. The tape search process looks at the first 
word in each block, therefore, the PLANETRACK program must con pute the CDN of 
the first word of the block containing the entry for the required day. 

Using as inputs the year, day of year, and the effect of integral and half-integral 
values of CELTIME, the program computes the Julian Day number for the typewriter. 
The ephemerides start and run continuously from 25 April 1963, which has a Julian 
Day Number of 2438144. 5; this number is subtracted from the computed Julian Day 
number and appears on tape as the Conventional Day Number. (it is 6 for | May 1963, 
the first day for which ephemerides tapes were made for the system. ) 

The first CDN in each block will be given by 6+ 32n where n is a positive 
integer or zero. If the first CDN of each record does not belong to this series, the 
tape search will be unsuccessful. ‘The first CDN in a block is stored in TAPEBLOCK 
(lower half) and the planet identification number is added to generate a tape search 
comparison word, stored in TAPESEARCH and also in SAFE. 

The number in SLOTBLOCK represents the number of the entry in the 32-word 
block. 

If 2 <SLOTBLOCK #=28, only one block need be read in; if not, either the pre- 
ceding or following tape block should also be read to provide sufficient ephemerides 
for 4th difference Bessel interpolation for a four-day period. 

If the one-or two-block search and read operations are successful, control is 
regained at NORMAL and the tape is rewound without interrupt or interlock, At this 
time, planet ephemerides for either 32 or 64 days are in core memory. 

The program performs 4th difference interpolation for Mercury and Venus and 
third difference interpolation for the other planets. 

Besselian interpolation of right ascension is done by the INTERPOL subroutine 
and the interpolated result is converted to hours, minutes, and seconds and stored 
for type out. Declination and distance are then interpolated and stored for type out. 
The other quantities which are to be typed out by INTERCOM are set up. Control is 
transferred to INTERCOM 7 times for the seven line type out, wnich completes oper- 


ation of the initialization section of PLANETRACK. 


B. Working Section 

The PLANETRACK working section begins at PLANETCONT. Control is trans- 
ferred to the subroutine DAY FIND which selects the day's ephemeris entry, stores 
the addresses of the table entries, and computes the interpolation argument P from 
the given value of CELTIME. The INTERPOL subroutine is then entered three times 
for the interpolation of right ascension, declination, radius and their time rates. 
An amount corresponding to the radius of the planet is subtracted from the interpo- 
lated value for accurate radar range control. All values are converted, scaled and 
stored properly in core memory and control is transferred to the master control 
program. 

C. Subroutines 

The initialization section of the PLANETRACK program makes use of the INTER- 
POL, DAYFIND, STATUSCK, LEFRNDOFF, and ROUNDOFF subroutines. The 
working section uses all except STATUSCK. 

Ll. INTERPOL 

The PLANETRACK interpolation subroutine uses Bessel's interpolation formiula 
(Reference 3) for 3rd or 4th difference interpolation of the ephemeris. The formula 
is 
ee ten + 6°f)] 


ft = I, - eer 4 (2 


=, Peath(Ps2) 1, ge eS) atod (P-2) (ot +0°t,) 


The quantities in the formula are associated with computer registers as follows: 
a. The address of f.. is in index register 4. 
b. The interpolation argument P is in register P. 
Cc Of) 19 is in register GAMMA. 
d. P(P-1) is in register PSQMP. 
e ot. - of, is in register DMINB. 


Bed) is in register PSQMP®O. 
3 


S 6 M2 is in register DMIN2CPLB. 


(o°8, zi 6") is in register EP2BM2DMA. 


rh 


aa 


ib. I is in the A register at the finish. 


The subroutine also performs numerical differentiation using the following formula 


(Reference 3): 


— 2P-1 
hf'p = Sf) .9 +— 


1 3p" Roar ae lea 6° 
————— 


2 2 
[ot OE 
172 


At the completion of the subroutine, hf’, is stored in the register NUMDERIV. 


When the subroutine is entered, Index Tascai 4 contains the address of fy 
and the register SETINTAD contains the address of ot of the quantity to be inter- 
polated. Index register 3 is used to acquire the other two first difference registers. 

2. DAYFIND 

The DAYFIND subroutine computes the Julian Day number and the addresses of 
the ephemeris entries corresponding to the current day. It also computes the interpo- 
lation argument P by properly scaling the sum of CELTIME and DELTATEE. The 
value of P is always positive, even though time be set negative by the system con- 
trol programs. The number in DAYINDEX, which can be -2, -1, 0, or +1, takes 
account of the effect of the integral part of CELTIME on the selection of ephemeris 
entry. 

3. STATUSCK 

The STATUSCK subroutine is entered from the external interrupt register asso- 
ciated with the tape channel when the interrupt occurs. Examination of the status 
code generates four possible outcomes: 

a. If the code indicates a normal completion, control is returned to 
the program at the interrupted point. 

b. If the code indicates that an end of tape mark was sensed, the tape 
is rewound and a message indicating failure to find the ephemeris is printed by the 
typewriter. The message is described in Section IIIB. Control is passed to the error 
return of the master control program. 

c. If the code indicates that an end of file mark was sensed, the regis- 


ter IMPERIAL is indexed and tested for the value 8. When less than 8 the next file 


is searched by passing control to TRYAGAIN. When IMPERIAL equals 8, an end of 
tape condition is assumed and the action is as described in (b) above. 

d. If the code indicates anything but end of file, end of tape, or normal 
completion, a tape error has occurred. The message described in Section IIIB, indi- 
cating the type of error is printed by the typewriter and control is passed to the error 
return of the master control program. 

4, LEFRNDOFF 
The LEFRNDOFF subroutine left shifts the AQ register the number of places 
indicated by the contents of index register 5 and rounds the A register. 
3S. ROUNDOFF 
The ROUNDOFF subroutine right shifts the AQ register the number of places 


indicated by the contents of index register 5 and rounds the A register. 
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